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H 1 ABSTRACT 

Context. Sequential supernova explosions create supershells which can break out a stratified medium producing bipolar outflows. The 
gas of the supershells can fragment into clouds which eventually fall toward the disk producing the so called galactic fountains. 
Aims. The aim of this paper is to calculate the expansion law and chemical enrichment of a supershell powered by the energetic feed- 
back of a typical Galactic OB association at various galactocentric radii. We study then the orbits of the fragments created when the 
supershell breaks out and we compare their kinetic and chemical properties with the available observations of high - and intermediate 
- velocity clouds. 

Methods. We use the Kompaneets (1960) approximation for the evolution of the superbubble driven by sequential supernova ex- 
plosions and we compute the abundances of oxygen and iron residing in the thin cold supershell. We assume that supershells are 
fragmented by means of Rayleigh-Taylor instabilities and we follow the orbit of the clouds either ballistically or by means of a hybrid 
5_( . model considering viscous interaction between the clouds and the extra-planar gas. 

Results. Given the self-similarity of the Kompaneets solutions, clouds are always formed ~448 pc above the plane. If the initial 
metallicity is solar, the pollution from dying stars of the OB association has a negligible effect on the chemical composition of the 
clouds. The maximum height reached by the clouds above the plane seldom exceeds 2 kpc and when averaging over different throwing 
angles, the landing coordinate differs from the throwing coordinate ~ 1 kpc at most. 

Conclusions. The range of heights and [O/Fe] ratios spun by our clouds suggest us that the high velocity clouds cannot have a 
Galactic origin, whereas intermediate velocity clouds have kinematic properties similar to our modeled clouds but overabundance 
observed for the [O/Fe] ratios which can be reproduced only with initial metallicities which are too low compared for those of the 

CO . Galaxy disk. 
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qq ' 1. Introduction related to the merging and accretion of these satellites by the 

f^i , Milky Way. In this picture, HVCs would be the gaseous counter- 

.. , Although significant progress has been made in the last few parts of the Milky Way circumgalactic stellar streams, which are 
< £ ■ years in exploring the distribution and chemical composition of be ii eV ed to represent the relics of dwarf galaxies that have been 
k> ; intermediate and high-velocity clouds (IVCs and HVCs, respec- accret ed by the Milky Way (e.g., Ibata 1994). While all these 
■ lively) in the halo of the Milky Way, an overall unified model models assume that HVCs ^ truly ex tragalactic objects which 
Ctf ■ for their formation is still lacking. These gas clouds are concen- m about to merge with the Galaxy from outside , there are other 
trations of neutral hydrogen (HI) with radial velocities which scenarios that see the IVCs and HVCs as objects that have their 
are not consistent with a simple model of differential galactic origin in the disk of the Milky Way , e.g., as pail of the so-called 
rotation (Richter et al. 2001). The distinction between HVCs » galactic fountains". In the galactic fountain model (Shapiro & 
and IVCs is loosely based on the observed radial velocities of Field 1976; Houck & Bregman 1990), hot gas is ejected out of 
the clouds; IVCs have radial velocities with respect to the Local the Ga i actic disk by supernova (SN) explosions, and part of this 
Standard of Rest (LSR) of 30 km s < | V LSS | < 90 km s while gas falls back in the form of condenS ed neutral clouds which 
HVCs have typical velocities | V LSS | > 90 km s . Several studies move at intermediate and high radial velocities. Whatever the 
(Lu et al. 1998; Wakker et al. 1999; Murphy et al. 2000; Bluhm origin of the Milky Way IVCs and HVCs is, it has become clear 
et al. 2001; Sembach et al. 2001) reveal different chemical com- that they must play an impor t ant ro i e i n the evolution of our 
positions for several of these clouds in different directions in the Galaxy. Two extremely important parameters to distinguish be- 
sky. Oort ( 1 970) proposed that these clouds represent condensed tween the Ga i act i c and extragalactic models of IVCs and HVCs 
gaseous relics from the formation phase of the Milky Way. This m the distance and me chemical composition of these clouds. In 
idea was later revived by Blitz et al. (1999), who suggested that case f the Galactic model the clouds ejected by SN explosions 
HVCs represent the building blocks of galaxies in a hierarchi- must have a velocity high enough to explain the observed dis- 
cal galaxy formation scenario. Since the Galaxy is surrounded tances and their chemical composition must be correlated with 
by smaller satellite galaxies (e.g., the Magellanic Clouds), an- the abundances of the gas in the disk 

other explanation is that IVCs and HVCs are gaseous streams In m pap£r we examine the scenario in which the Qrigin 
of these clouds is given by SN explosions, as the Galactic foun- 
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physical and chemical properties of the swept up gas by sequen- 
tial SN events, in literature called supershell. Once the clouds 
leave the stellar disk, they move through the extra-planar gas 
halo. Two extreme types of models have been considered here 
for the motion of the gas through the extra-planar gas: i) the 
ballistic and ii) the fluid homogeneous models. Numerous non- 
homogenous filamentary structures observed in the extra-planar 
halos (NGC891, NGC 5775) suggest that the kinematics of these 
halos could be interpreted with ballistic models. Ballistic mod- 
els describe the gas as an inhomogeneous collection of clouds, 
subject only to the gravitational potential of the galaxy: for ex- 
ample, in the galactic fountain model the ejected gas from SN 
events falls back ballistically (Bregman 1980). These models are 
able to explain vertical motions of the cold and warm gas com- 
ponents observed in several spiral galaxies (e.g. Fraternali et al. 
2004; Boomsma et al. 2005). However, Collins et al. (2002) have 
tried a ballistic model to describe the extra-planar ionized gas of 
NGC 891 and found problems in reproducing the observed kine- 
matics. The ballistic model predicts that clouds migrate radially 
outward as they cycle through the halo. The mass fluxes esti- 
mated from the models of NGC 891 and NGC 5775 imply that 
significant amounts of gas can be involved in these migrations. 
Such migrations could cause a redistribution of gas that could 
affect metallicity gradients as well as star formation properties. 
Such effects have been previously investigated by Charlton & 
Salpeter (1989), for example, but extensive observations of the 
kinematic behavior of edge-on galaxies should yield important 
constraints on such redistribution. Bregman (1980) presented a 
model where the clouds ejected fall into the disk subject to the 
gravity of the Galaxy at the same coordinate where they were 
thrown. In this case the metallicity gradient of the Galaxy does 
not change. 

Collins et al. (2002) suggested that the ballistic model prob- 
lems could be solved by considering the presence of drag be- 
tween disk and halo, such as through magnetic tension or vis- 
cous interactions between clouds. Barnabe et al. (2006) pre- 
sented fluid stationary models able to reproduce the observed 
negative vertical gradient of the rotation velocity of the extra- 
planar gas. In these models the gas is described as a stationary 
rotational fluid. In the case of a single cloud, the solutions of this 
model predict the drag invoked by Collins et al. (2002), as ram 
pressure between cloud and fluid halo. This suggests that a cor- 
rect description of the extra-planar gas dynamics may be found 
in hybrid ballistic-fluid stationary models. Recently, Fraternali 
& Binney (2008) presented a new formulation for viscous cloud- 
halo interactions in the framework of the cold gas accretion onto 
spiral galaxies in the local universe. 

In this paper we will follow the clouds ejected by SN explo- 
sions both in the Galaxy in the purely ballistic model and in the 
hybrid one. The paper is organized as follows: in Sect. [2] we will 
describe our model, in Sect|3]we will show our results concern- 
ing the chemical abundances of the supershells and the cloud 
orbits in the purely ballistic and in the hybrid model. In Sect. [4] 
we compare our results with observations of the HVC Complex 
C and the IVC VIArc. In Sect. |5]we report our conclusions. 

2. Description of the model 

Type II SNe usually occur in OB associations containing sev- 
eral dozens of massive stars. Sequential SN explosions create 
a superbubble, whereas the swept up gas is concentrated in a 
thin cold shell called supershell. The superbubble expansion in 
a stratified medium does not follow a spherical evolution and 
the Kompaneets (1960) approximation well describes the way in 



which the superbubble grows up in the meridional Galaxy plane. 
After some interval of time, a region of the supershell can frag- 
ment due to the occurrence of the Rayleigh-Taylor instabilities 
(RTIs), therefore clouds of gas can form. Once left the stellar 
disk, the orbit of each cloud can be followed either ballistically 
or with a hybrid model considering viscous interaction between 
the cloud and the extra-planar gas. 

2.1. The Kompaneets approximation 

The physical parameters of the superbubble can be fairly well 
described by the Kompaneets (1960) model. The Kompaneets 
approximation assumes the following : uniform pressure within 
the superbubble, superbubble expansion in a direction normal 
to the local surface and expansion speed implied by a strong 
shock (i.e., the internal pressure dominates the external pres- 
sure). Kompaneets found an analytic expression for the shape 
of the bubble during its expansion in an exponential atmosphere 
with density: 



p(z) = p exp(-z/H), 



(1) 



where po and H are the disk density and height scale, respec- 
tively. We consider the shock expansion in cylindrical coordi- 
nates (z, R). In the work of Kompaneets it was assumed that for 
z < the ISM was denser than the gas residing above the disk. 
Therefore, in this approximation, the bubble has an asymmetric 
evolution with respect to the plane of the galaxy. Kompaneets 
found analytic expressions for the top and bottom sites of the 
remnant. Referring to Fig. Q] we have: 



z UH -2mn{i- + 2.), 

and 



h - 2/7 arcsm (^ : ) ■ 



(2) 



(3) 



The quantity y is a transformed variable (with units of length) 
defined by: 



v 



Jo \ 2 poQ 



(4) 



where Q. and £,/, and y are the volume, the thermal energy and 
the ratio of specific heats for the superbubble, respectively . The 
thermal energy is calculated from the differential equation: 



dE t h 
dt dt ' 



(5) 



where Lq is the wind luminosity and the pressure is given by the 
equation of state: 



(6) 



In our work for the ISM z-profile we assume an exponential law 
also for z < 0. Hence the superbubble evolves symmetrically 
with respect to the galactic plane resulting in a peanut-like shape, 
as suggested by Mac Low et al. (1988) and by Tenorio-Tagle et 
al. (1999). The volume of the superbubble can be estimated by: 



Q = 2 • ^na 2 b - ^n\z H \ 3 



(7) 



where a - {zl + \zh\)/2 is the semi-major axis. A numerical in- 
tegration of eqs. (H]i and ©, with the help of eqs. (|6]l and (0, 
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Fig. 1. Explanatory picture of the superbubble shape (see 
Kompaneets 1960) under the assumption that the density below 
the plane is constant (see the Text). There are reported some 
physical quantities: the height of the superbubble in the expo- 
nential atmosphere (Z L ) and below tha plane (Z H ), and the semi- 
minor axis (b). Figure taken from Kompaneets (1960). 



yields y(t) and E t h which implicitly give the time evolution of 
the superbubble. The solution depends on only three parameters: 
the scale height H, the density near the source po, and the wind 
luminosity Lq. Hence the problem is solved most naturally in a 
system of units which are determined by these three parameters. 
The unit of length is then Hq, the unit of mass is poH 1 , and the 
unit of time is (poH 5 /Lo) l/3 . We obtain the dimensionless solu- 
tion by integrating the dimensionless form of eqs. (O and (0. 
In Figs. [2] and [3] we report the dimensionless evolution of y and 
the velocity dzL/dt. The dimensionless analysis predicts that the 
superbubble begins to accelerate at (see Fig. [3): 



bitrary but in this way our results are consistent with the work of 
Mac Low & McCray (1988). We have also tested velocity varia- 
tions of 5% and 20% and the results are not very different to the 
ones presented in this paper. In our model we assume that ow- 
ing to the RTIs, the supershell fragments and we consider each 
fragment as a cloud with initial velocity given by the supershell 
velocity at the moment of the fragmentation. RTIs arise when the 
supershell accelerates since at that moment, the effective gravity 
of the supershell exceeds the one of the tenuous superbubble. 
From Tenorio-Tagle et al. (1987) we have taken an estimate of 
the time of growth of these highly irregular structures: 



tR-T 



io 5 /4 



no 



£51 (fa + ?Ar) 



no yr, 



(9) 



where A\q and R40 are perturbation wavelength and remnant ra- 
dius measured in units of 10 and 40 pc, respectively, £51 is the 
energy measured in units of 10 51 erg and t& v is the time necessary 
for reaching an appreciable acceleration. At the time when RTIs 
start growing, the supershell has reached z/,=3 12 and Z?=228. We 
assume for simplicity the formation of a single cloud that cor- 
responds to the maximum perturbation wavelength for the thin 
shell, roughly given by the supershell thickness. In the case of 
adiabatic expansion of a spherical thin shell, the thickness 6 is 
given by 5 — R/12, where R is the radius of the shell. We as- 
sume for simplicity that 6 is the maximum perturbation wave- 
length and we make the approximation that R = zl, therefore, 
6 — A = 26 pc. If we define Ir-t as: 



tR-T - 



tR-T 



(p // 5 /Xo) 1/3 



(10) 



we obtain the total adimensional time necessary for the growing 
of instabilities and for the fragmentation of the bubble: 



tfinal = t a + ?Av + tR-T ~ 4.37, 



(ID 



Considering different values of the parameters p , H, Lq, we ob- 
tain different times at which the cloud could leave the disk: 



tfinal =4.37(p H 5 /£o) 



1/3 



(12) 



Since the internal pressure dominates the external pressure, the 
expansion speed is that given by the Hugoniot condition for a 
strong shock: 



7+1 POO 



\ 2 p{zY 



(13) 



ta = 



(p„// 5 /£o) 1/3 



1.72, 



(8) 2.2. Abundances of Fe and O in the superbubble 



namely when dzi/dt has a minimum. We underline that the sim- 
ulated superbubbles, characterized by different values of po, H, 
Lo, reach the accelerated phase at different times, as shown in 
eq. (O, but they have the same shape and volume. This is a di- 
rect consequence of the self-similar nature of the Kompaneets 
solutions. Using the scale height H-141 pc (see Sect. 12.3b we 
obtain for the height of the supeshell and for the minor semiaxis 
the values zl = 191 pc and b = 145pc. From the Fig.[3]we note 
that, after the beginning of the accelerated phase, the velocity 
remains roughly constant for ~ 2 t. In this interval of time the 
acceleration is therefore negligible. Therefore we assume that 
the RTIs begin to grow when the velocity is increased by a fac- 
tor of 10%, because from that moment on we can assume that the 
acceleration is not negligible. The choice of this percentage is ar- 



In this section we compute the total amounts of O and Fe that 
reside in the supershell at the time f/,„ fl / (see eqlTZb in which we 
suppose the fragmented shell can leave the stellar disk of the 
Galaxy. We assume that all the stars with masses larger than 8 
M explode as Type II SNe and we adopt the following main 
sequence lifetimes (Matteucci & Greggio 1986): 



1.2m" 1 ' 85 + 0.003Gyr. 



(14) 



In our model, the upper mass limit of the OB association is as- 
sumed to be 4OM . The lifetime for a 4OM is, from eq. ([T4i >. 
t 40 - 4.30 x 10 6 yr. 

For any parameter choice, at the time in which the shell frag- 
ments and leaves the disk, the smaller mass stars in the OB as- 
sociation are still alive. Using eq. (TBI we are able to estimate 
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Fig. 2. Evolution in the dimensionless space of the quantity y 
(eq. H]i in function of t . y is related to the top of the superbubble 
(Zl) evolution from the eq. @. The dimensional time t is given 
by: t = t x (p # 5 /£o) 1/3 



Fig. 3. Evolution in the dimensionless space of the velocity 
dfi/dt as function of t. We note that the shell starts accelerat- 
ing at t a — 1 .72 . 



the lower mass limit for a SN event corresponding to the time at 
which the cloud forms. The time at which the cloud is formed is 
given by ?/,■„„/ plus the lifetime of the largest considered mass, 
namely: 



Tlim = t final + T"4Q. 



(15) 



The mass of metals (in our model we consider only O and Fe) 
ejected by the range of SNe considered is given by: 



/~40 



m e i(m)(f>(m)dm, 



(16) 



where <p(m) is the Initial Mass function (IMF) assumed to be the 
Salpeter (1955) one, namely: <p{m) = Am~ (2 35) . A is the normal- 
ization constant, given by: 



I <p(m)dm — A 



m' 235 dm = SNe, 



(17) 



where SNe is the assumed number of SNe in the OB association. 
In eq. ( fT6] l the term m e i(m) is the ejected mass of the el element 
considered (O or Fe) by the star of initial mass m (Woosley & 
Weaver 1995). We consider the stellar yields at four different 
metallicities: solar metallicity Z , 0.1 Z , 0.01 Z , 1O~ 4 Z . We 
assume that the ISM in the disk has the same metallicity of the 
OB association. For the solar abundances (by mass) we use the 
Anders & Grevesse (1989) values: 



Oi 6b = 9.59 x 10- 3 , Fe 56o = 1.17 x 1(T 



(18) 



We use this set of solar abundances in order to be consistent with 
the yields of Woosley & Weaver (1995) which refer to Anders & 
Grevesse (1989) abundances. However it should be recalled that 
recently Asplund et al. (2005) recomputed the O solar abundance 
and deduced a log(O/Fe) lower by ~ 0.24 dex than the Anders 
& Grevesse (1989) one. The total amount M s heii-el of the mass 
of the element which resides in the thin shell of the superbubble 
is given by the eq.: 



M 



shell- el 



M 



shell ■ 



Z + M, 



(19) 



where Z is the initial metallicity of the OB association (and con- 
sequently of the disk ISM). We adopt the extreme approximation 
that total amount of the ejected metals from SNe ends up in the 
shell. M S heii is me mass of the ISM swept up into the thin shell 
and is given by: 



M she „{z > 0) = np b 2 



-m 



z-a + \z H \ \ 



2\ 



dz. (20) 



We underline that for the Fe we use only solar yields for all 
range of metallicities (Chiappini et al. 2001). 

2.3. Galaxy model 

The potential well of the Galaxy is assumed to be the sum of 
three components: a dark matter halo, a bulge and a disk. The 
dark matter halo gravitational potential is assumed to follow the 
Navarro, Frenk and White (1996) profile: 



-4nr 2 



'dmfi Pdmfl " 



ln(l +x) 



(21) 



where pd m ,o is a reference density, r^ m o is a scale radius, x = 
r/rdm,o and r is the spherical radius. The halo is truncated at a 
radius r^mt beyond which its potential follows the 1/r profile. 
The bulge gravitational potential is given by (Hernquist 1990): 



®b(r) = - 



GM b 
rb,o + r 



(22) 
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where r^o is a scale radius and Mb is the bulge mass. For the 
disk potential we have chosen the axisymmetrical Miyamoto & 
Nagai (1975) model, which in cylindrical coordinates (R,z) can 
be written as: 



GM d 



^R 2 + (a + Vz 2 + b 2 f 



(23) 



Tab. (0 gives the values of all the parameters concerning the 
Galaxy model. For the Interstellar Medium (ISM) z-density pro- 
file we used the eq. (HJ where: po = nopm p is the density in 
the disk plane; m p is the proton mass and /i is mean molecu- 
lar weight for the disk (assumed to be 0.61). For the value of 
H we adopted the vertical distributions of the various interstel- 
lar components (molecular, cold HI, warm HIa, warm HIb, HII 
Regions, and diffuse HII) in the solar neighborhood as reported 
by Cox (2005). We obtained: 



H= — V Pi {z)dz * 141pc. 
Po Jo j-* 



(24) 



The ISM density profile along the radius R of the Galaxy is taken 
by Wolfire et al. (2003). In our model we consider the presence 
of an extra-planar gas halo. This halo can be described by means 
of the perfect gas law: 



pkT 
yan v ' 



(25) 



For simplicity we assume a isothermal and self-gravitant distri- 
bution, hence introducing the parameter: 



kT 

A> = — , 

fim p 



(26) 



by means of which the eq. d25l ) can be written as: 

P=PoP. (27) 

Using the hydrostatic equilibrium equation with the assumption 
of a static halo of gas we have: 



p(R, z) = p e' 



(-av+4>o)/A) 



(28) 



where <P, , = <S>d m + <D>£ + <E>m and Oo is the value of $> tot calcu- 
lated in the galactic center. The temperature T a has been chosen 
following the relation: 



\kT B 



3m p pv 2 c 2 



(29) 



where v C 8 is the circular velocity in the plane of Galaxy in the 
solar neighborhood and we obtain that: 

T = 3.7 x 10 6 K (30) 
for the chosen parameters. 

2.4. Galactic fountains in a ballistic model 

Once the top of the supershell reaches the height above the galac- 
tic plane related to the time ?/,■„„/ (see eq.[T2]i. the thin shell can 
leave the stellar disk and move towards the extra-planar gas. 
Ballistic models describe the gas as an inhomogeneous collec- 
tion of clouds, subject only to the gravitational potential of the 
Galaxy. The fragments of the shell have different initial velocity 



modulus || v ||= v„ and masses depending on the the number 
of SNe in the OB association and the initial throwing radial co- 
ordinate. We have chosen for our simulation Cartesian coordi- 
nates with versor (e x , e v , e,-). Given a generic point in this space 
(x, y, z), it corresponds in the meridional plane of the Galaxy to 
(R, z) — ( yjx 2 + y 2 , z), where R is the radial coordinate. Since the 
|| \ || velocities are relative to the local standard of rest for the 
trowing radial coordinate Rq, we have that in the inertial refer- 
ence frame of the simulation in Cartesian coordinates: 



Vinitiai(Ro) = v a + v c (R , 0)e y . 



(31) 



For each choice of || \ II we consider several throwing direc- 
tions to take into account the possibility that the breakup of a 
single supershell produces several fragments, each one with its 
own orbit. Given a throwing direction, a generic vector Vo has 
components along e x , e z e e v with respect to our inertial frame: 



v zm — v„cosy 
Vxj„t, - v„sinycos/? 
Vy M = v„sinysinyS + v c (Rq, 0) . 



(32) 



The parameters y e (3 vary in order to recreate, in the local stan- 
dard of rest, a symmetrical fountain of jets: 



Pi = — with; = 0, 1...7 



and 



^2-177 



with ,/ = 2,3, ..5. 



(33) 



(34) 



We consider also the case in which vo has only a component 
along e~. In conclusion, for each velocity modulus v„ we calcu- 
late a fountain composed by 33 jets. 

2.5. Beyond the ballistic model: hybrid ballistic-fluid 
stationary model 

In this paper we want to provide some hydrodynamical modifi- 
cations to previous ballistic model. In the motion equations for a 
single cloud we insert a Stokes term in order to consider also the 
viscous interactions between cloud and extra-planar gas. We in- 
troduce the drag time td rag , i.e the time required for an HI cloud 
to corotate with a homogeneous gas halo (Barnabe et al. 2005): 



3Cd Vrd 



(35) 



where Co — 1 is a numerical coefficient, r c is the radius of a typ- 
ical H I cloud, v r ei is the modulus of the relative velocity between 
the cloud and the homogeneous extra-planar gas, and^ = p c /p 
is the ratio between the cloud and the medium densities. td rag 
was estimated assuming pressure equilibrium between cold and 
hot components, from which x — T/T c m 3000. A fiducial value 
of the relative velocity is obtained by assuming v re i ~ 2 x 10 7 
cm while the cloud radius r c is estimated as: 



4nr c 3 pJ3 = 1O 5 M M 5 , 



(36) 



where Ms is the cloud mass in units of 10 5 M Q , thus following 
the work of Barnabe et al. (2005) one obtains: 



rag 



2.7 x 10 s 



1/3 



(37) 



6 

Table 1. Galactic parameters. 



Spitoni et al.: Galactic fountains 



Pdm.O 


M b 


M d 


rdm.O 


f. dm J 


''/>.(> 


a 


b 


1(T 24 g cm - 3 


10 10 M o 


10 10 M o 


kpc 


kpc 


kpc 


kpc 


kpc 


0.29 


3.5 


7.69 


30.8 


347.7 


0.8 


8.45 


0.26 



where n p is the numerical density of the extra-planar gas with Table 3. Same of Tab.[2]but for Rq = 8 kpc. 
the density profile described by eq. ((28). We note that a viscosity 
term could not be directly inserted in a Lagrangian formulation, 
thus we must write the explicit drag term to add in the motion 
equations. The drag term can be derived from the brake force 
estimate as: 



a dras = (38) 

''drag 



SNe 


* final 


v„ 




[yr] 


[kms" 1 ] 


10 


1.90x10' 


27 


50 


l.llxlO 7 


46 


100 


8.84xl0 6 


58 


500 


5.17xl0 6 


100 



where v is the velocity of the cloud whereas v g is the velocity for 

the extra-planar gas that we assumed equal to zero. Table 4. Same of Tab.[2]but for R = 12 kpc. 



3. Our results 

3. 1 . Chemical composition of the clouds 

In our models we vary the number of SNell in the OB associ- 
ation (SNe) and the throwing radial coordinate (Rq). We con- 
sider four possible OB associations containing 10, 50, 100, 500 
SNe respectively. Assuming an explosion energy of 10 51 erg, 
the luminosities Lq of these OB associations are 10 37 , 5 x 10 37 , 
10 38 and 5 x 10 38 ergs -1 respectively. These numbers of mas- 
sive star in OB associations are consistent with the observations 
(de Zeeuw et al. 1999). We simulate galactic fountains with 3 
different throwing radial coordinates Rq: 4, 8, 12 kpc. Given the 
assumed Galaxy model, Rq defines the disk density po, whereas 
the scale height is constant (see Sect I2.3l l. In Tabs. [2] [3] [4] we 
summarize the results for fragmentation times and the velocities 
of the superbubbles in the direction perpendicular to the galactic 
plane at those times for 4 kpc, 8 kpc, 12 kpc, respectively. 

Table 2. Cloud formation times and cloud velocities for Rq = 4 
kpc. 



SNe 


* final 


v„ 




[yr] 


[kms" 1 ] 


10 


2.29x10' 


23 


50 


1.34X10 7 


39 


100 


1.06xl0 7 


49 


500 


6.21xl0 6 


83 



For all the values of SNe and Rq considered at the time 
t fi na i, at which clouds are thrown out of the disk, the supershell 
presents: zl = 448 pc, zh = 165 pc and b = 259 pc. Our results 
are in agreement with the work of Mac Low & McCray (1988): 
instabilities, for roughly the same luminosity range, become im- 
portant at 3H height scale. The total mass of gas swept up by 
the SN shock wave for positive z-coordinates is given by the eq. 

C39. 



SNe 


* final 


v„ 




[yr] 


[kms- 1 ] 


10 


1.67x10' 


31 


50 


9.77x10 s 


53 


100 


7.74x10 s 


67 


500 


4.53x10 s 


114 



In conclusion, we obtain that the masses of the ISM swept 
up into the thin shell for Rq= 4, 8, 12 kpc, respectively are: 



M 4 (z > 0) = 10.07 x 10 s M 0: 
M 8 (z > 0) = 5.79 x 10 5 M o , 
M n (z > 0) = 3.89 x 1O 5 M . 



(39) 
(40) 
(41) 



All the results about the O and Fe abundances in the clouds 
ejected by sequential SN explosions as functions of SNe and 
Ro are reported in the Online material. M,f e5b and M«o 16 are the 
total amounts of Fe^e, and Ok, in unit of M , whereas M+ e j is 
the total mass ejected by the OB association: 



M 



*ej 



JM inf 



m, ot (m)<p(m)dm, 



(42) 



where m lol {m) is the total mass ejected by a SN as a function 
of its initial mass and metallicity. M tot , X* Fe56 , X, 0m and [O/Fe] 
are: 



M, ot = M sheU + M+ej 
MshellFe 56 



X, 



M„ 



. x *o l6 - 



M. 



SheWht 
Mm 



[O/Fe] = log 



MshellO l6 
MshellFe 56 



log 



O 



16, 



Fe 56e 



(43) 
(44) 

(45) 



In Fig|4]we show the predicted [O/Fe] ratio as a function of the 
number of SNe and of the initial metallicity in the solar vicinity. 
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In the meridional plane of the Galaxy the initial conditions are 
(R, z)=( 8 kpc, 448 pc). We note that significant over-abundances 
of O relative to Fe are found only in the case of a large num- 
ber of SNe and low initial metallicity. In Fig. [5] we report the 
same quantities but using stellar yields given by Kobayashi et 
al. (2006). In Fig. [6] we show [O/Fe] ratios as functions of the 
number of SNe assuming solar metallicity but varying the initial 
throwing coordinate R a . We note that larger radial coordinates 
yield a larger [O/Fe], because the amount of the swept-up pris- 
tine gas is smaller (see eq.l4TT> and therefore the new a elements 
ejected by SNe are less diluted. In Fig. [7] we report the [O/Fe] 
ratios varying the initial throwing coordinate Rq and taking for 
the initial ISM metallicities the average observed values given 
by Andreievsky et al. (2002a-c, 2004) and Luck et al. (2003), 
as a function of galactocentric distance, by analyzing Galactic 
Cepheids (see Cescutti at al. 2006). Referring to the Tab. 4 of 
the work of Cescutti et al. (2006), we find: Z = 1 .65 x Z for 
R = 4 kpc and Z = 0.74 x Z Q for R = 12 kpc. 




N SNe 



Fig. 4. [O/Fe] ratios as functions of the number of SNe and of 
different disk gas metallicities in the case of Ro=8 kpc using 
stellar yields given by Woosley & Weaver (1995). 



3.2. Dynamics of the galactic fountains 

Our analysis focuses first on the study of solar neighborhood 
galactic fountains (i.e. with Rq = 8 kpc). In Fig. [8] we show our 
results concerning the orbits of the shell fragments once they 
leave the stellar disk in the purely ballistic model. The trajec- 
tories of the orbits of each galactic fountain are plotted in the 
meridional plane of the Galaxy {R,z). Spatial initial conditions 
are the same for all cases considered: Rq = 8 kpc and zo = 448 
pc, whereas the initial velocities depend on the number of SNe 
considered (Tab. [3). In Fig. [8] squares on the R axis reported the 
average falling radial coordinate in the Galaxy plane. As we can 
see, the clouds are preferentially thrown outwards, but their fi- 
nal average landing coordinates differ by 1 kpc at most from the 
throwing coordinate. This result is consistent with the works of 
Bregman (1980), Fraternali & Binney (2008) and Melioli et al. 



Kobayashi el al.(2006) 


■ 






0.2Z e 




0.05Z o 

— M 

_____ . 




10 50 100 

N SNe 


500 



Fig. 5. As in Fig|4]using stellar yields given by Kobayashi et al. 
(2006). 



0.03 




10 50 100 500 



N SNe 

Fig. 6. [O/Fe] ratios as functions of the number of SNe and the 
throwing radial coordinate in the case of solar metallicity using 
stellar yields given by Woosley & Weaver (1995). 



(2008). On the other hand, in their recent hydro-simulations for 
the Milky Way disk Booth & Theuns (2007) found that clouds 
ejected from galactic fountains return to the disk at average 
galactocentric distances several times larger than the galacto- 
centric distance of the fountain. However, their mass resolution 
(particle mass 1.5 x 1O 5 M in the highest resolution simu- 
lation) makes their results more suitable for understanding the 
global behavior of the HI in the Galaxy rather than the evolution 
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Fig. 7. [O/Fe] ratios as functions of the number of SNe and the 
initial throwing radial coordinate, taking for the ISM metallic - 
ities the average observed values along the Galactic disk from 
Cepheids by Andreievsky et al. (2002a-c, 2004) and Luck et al. 
(2003) (see Cescutti et al. 2006). 



of a single superbubble/supershell. In Tab. |5]we report the main 
radial velocity for each fountains given by: 

< AR > 



< V R > = 



< At > 
where AR - R 



final 



- R and At - t 



final 



+ t, 



orbit 



the time f. 



(46) 



rbti IS 



the time required for the cloud to return to the galactic disk once 
it leaves the supershell. Several chemical evolution studies have 
suggested that metallicity gradient disappears if the velocity of 
the radial flow is |v/| > 2 km s" 1 (see Matteucci 2001). 

Table 5. Radial flow velocities (in km s~') in the solar neighbor- 
hood (model without drag) as a function of the radial throwing 
coordinate. 





4kpc 


8kpc 


12 kpc 


10 SNe 


0.6 


0.6 


1.1 


50 SNe 


1.6 


2.3 


4.2 


100 SNe 


2.6 


4.1 


6.6 


500 SNe 


4.5 


13.2 


14.4 



For the model with viscous interaction we must know also 
the mass of the cloud because the drag terms depend on this 
quantity (see eq.[37l). Referring to Fig.Q]we assume that the part 
of the shell that could fragment and skip upwards is the mass 
included above the s height where 4| = 0. Over this height each 
point of the shell has a positive velocity component along the 
z axis. Therefore, we estimate that the masses of a gas cloud 
ejected by sequential SN events in the extra-planar halo are of 
the order of: 



Fig. 8. Galactic fountains reported in the meridional plane in the 
purely ballistic model with the same spatial initial conditions: 
(R,z)=(8 kpc, 448 pc). Squares on the R axis are the average 
falling radial coordinate. 



M, 



doud « 8 



loudn =v. 



1.24 x 10 3 M o 



: 0.84 x 10 3 M R 



(48) 



(49) 



Mdoud^ =2.17x 10 5 M G 



(47) 



We note that the masses of the clouds are roughly 20% of the 
initial supershell mass, therefore the total amount of the ejected 
metals in clouds are roughly 20 % of the produced metals by 
the OB association, in agreement with the recent Melioli et al. 
(2008) hydrodynamical simulation results. 

In Fig.|9]we show that for the range of initial velocities con- 
sidered here, the effect of a viscous term in the motion equations 
is weak. The most evident effect is reported in Fig. [10] where 
the hybrid model is compared with the purely ballistic one in 
the case of 500 sequential SN explosions at Rq— 12 kpc. The 
natural effect of a viscous interaction between the cloud and the 
extra-planar gas is to brake the motion of the clouds. Therefore, 
the average radial falling coordinate on the disk for the model 
with drag < R/tnai >= 12.37 kpc is much smaller than the one 
obtained with the purely ballistic model (< R final >= 14.38 kpc) 

As shown in Tab. [6] it is likely that the most realistic number 
of massive stars in OB associations in our Galaxy is about 100. 
In Fig. QT|we have reproduced various fountains for the model 
without drag obtained with 100 sequential SN explosions respec- 
tively at 12, 8, 4 Kpc as radial radial coordinates. As shown be- 
fore, different radial initial conditions imply different ejection 
velocities for the clouds. In particular, the throwing coordinate 
Ro-12 kpc is surrounded by a more tenuous disk, therefore the 
velocity is larger and the orbits wider. The highest z coordinate 
that a cloud could reach in this case is about z= 2.3 kpc. 

4. Comparison with HVCs and IVCs observations 

In this section we compare our models with observational data 
concerning distance and chemical composition of some clouds in 
our Galaxy. We consider two cloud systems: Complex C (HVCs) 
and IV Arch (IVC). 
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Fig. 9. Galactic fountains in a model with drag with the initial 
radial coordinate Ro=8 kpc. In this case the orbits depend on the 
mass of the cloud ejected (indicated on the top of the panels). 
The natural effect of the viscous interaction is to decelerate the 
clouds. 



114 km g- 1 500 SNe M D = 0.83 xlO 5 M s 
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Fig. 10. Fountains which reach the maximum height in our 
model. 

4.1. Complex C 

In the work of Wakker et al. (2007) it was reported the first suc- 
cessful detection of interstellar Ca n H and K absorption from 
HVC complex C. They concluded that complex C is located at 
Galactocentric radius <14 kpc, and lies high above the Galactic 
plane (z=3-9 kpc). 

Integrating the Hi column density across the cloud, it was 
estimated M( Hi) as 0.7-6.0 xlO 6 M . As we show in eqs. 
d47i i. (US), (|49l , the sizes of our ejected clouds are considerably 
smaller. Even if we consider the case with the superbubble given 



Fig. 11. Fountains in the purely ballistic model driven by a se- 
quential explosion of 100 SNe 



Table 6. Massive stars in OB associations from P.T de Zeeuw et 
al. (1998) 



OB association 


WR 


O 


B 


A 


Total 


Upper Scorpius 






49 


34 


83 


Upper Centaurus Lupus 






66 


68 


134 


Lower Centaurus Crux 






42 


55 


97 


Vela OB2 


1 




81 


5 


87 


Trumpler 10 






22 


1 


23 


Collinder 121 


1 


1 


85 


8 


95 


Perseus OB 2 






17 


16 


33 


a Perseus OB 2 






33 


30 


63 


Lacerta OBI 




1 


35 


46 


82 


Cepheus OB2 




1 


56 


10 


67 



by 500 SN explosions with Rq-12 kpc, the highest z we can 
reach with this event is about z — 4.4 kpc and at this height 
our cloud velocities are roughly equal to zero. In the Tab. 1 of 
Wakker et al. (2007) are reported the velocities (relative to the 
LSR) of complex C in the direction of different line of sight, and 
in average the complex C velocity is -114 km/s. Therefore, the 
kinematical data of Complex C are inconsistent with our results. 

A key point for understanding the origin of this HVC is con- 
sidering its chemical composition. Because oxygen is not sig- 
nificantly depleted onto dust grains (Meyer et al. 1998) and the 
Oi/Hi ratio is not altered by ionization effects, we have that 
(Oi/Hi) * (O/H). In the work of Richter et al. (2001) the so- 
lar values for the abundances are given by Anders & Grevesse 
(1989) and Grevesse & Noels (1993). For the oxygen they ob- 
tained: 

[O/H] - [O i/Hi] = -L03!<>^. (50) 
For the iron abundance in the Complex C it was obtained: 
[Fe/H] * [(Fe n + Fern)/ Hi]. = -1.15, (51) 
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Based on the data from Fe n they obtained [Fe/H] - -1.27^2 j9. 
Therefore a low iron abundance in Complex C could be a result 
of depletion of iron into Fe rich dust (Savage & Sembach 1996). 
A more careful discussion about the possible effect of dust is in 
Sect. 4.3. Hence the relative abundance is [O/Fe] = 0.12. In our 
model a similar result is obtained in the case of 500 SNe with 
initial radial throwing coordinate fixed at Rq = 12 kpc and, as 
disk initial gas metallicity, Z = 0.1 x Z . As we have seen, the 
most likely metallicity of the disk at Ro= 12 kpc at the present 
time is 0.74 Z , which is not consistent with the total metallic- 
ity of Complex C (eqs. 50 and 51). We obtain with our model 
[O/Fe] = 0.03 dex in the case of 500 SNe with throwing co- 
ordinate = 12 kpc and metallicity fixed at 0.74 Z (see Fig. 
[7]). Such an oxygen overabundance, although not impossible, is 
quite unlikely. Given the inconsistency of the kinematical data 
with our predictions, we can rule out a Galactic origin for the 
Complex C HVC. 

4.2. IV ARCH 

Richter et al. (2001) suggested that IV Arch, given its nearly 
solar abundance and its z-height bracket of 0.8-1.5 kpc, has its 
origin in the Milky Way disk, probably as part of a Galactic 
fountain. These heights are indeed consistent with our results, 
and also the velocities of IVCs are easily understandable in the 
framework of galactic fountains. The abundances measured by 
Richter et al. (2001) are: [O/H] = -0.01^ and [Fe/H] = 
-0.26+° }9, therefore: 

[O/Fe] = [O/H] - [Fe/H] = 0.25 dex. (52) 

If we assume a Galactic fountain origin, our model yields 
[O/Fe] = 0.24 dex in the case of 500 SN explosions, R = 8 kpc 
and metallicity Z = 0.01 xZ , and [O/Fe] = 0.26 dex in the case 
of 100 SN explosions, Rq= 12 kpc and metallicity Z = 0.01 xZ . 
As we have seen, it is quite unlikely that the initial metallicity of 
an OB association can be nowadays as low as 0.01 Z and this 
low initial metallicity would be inconsistent with the nearly solar 
abundance measured in IV Arch. Therefore, although the kine- 
matics properties of IV Arch are in agreement with our results, 
its large [O/Fe] casts some doubt about the Galactic origin of 
this IVC. 

4.3. The effect of dust 

In the work of Lu et al (1998) it was shown how important is 
the depletion of metals (in particular Fe) in the ISM into dust. 
Since this can affect out results about [O/Fe] ratios, here we 
estimate the time scales of the destruction and accretion of dust 
in a superbubble. Referring to the work of Calura et al. (2008) 
we define, for a given element i the destruction time-scale Tdestr 
for the dust in the superbubble due to SN shocks as: 

Tdestrj = (eM SNR ) ■ — — , (53) 

1<SN 

where R$n is the total SN rate, M$nr is the mass of the inter- 
stellar gas swept up by the SN remnant and e is the destruc- 
tion efficiency in a three-phase medium. We consider for ctism, 
the surface gas density of the ISM, the value of 1OM pc~ 2 and 
cMsnr — 1O 5 M . For the SN rate we consider the typical case 
of the OB association with 100 SNe. We assume as radius of 
the OB association the conservative value of 100 pc. In the end 
we obtain that Tjestr.i — 0-7 Myr. This timescale is considerably 



smaller than the time necessary for the formation of a RT in- 
stable supershell in our model (See sect. 13. II ). If the number of 
SN is smaller, Tdestr is larger but it remains considerably smaller 
than t fi na i. Given the complexity of the interaction between SN 
shocks and dust (Draine & Salpeter 1979; Jones et al. 1996) this 
estimate is necessarily simplified but it leads us to the suggestion 
that only a small amount of dust can survive in the supershell 
swept up by the OB association, therefore only a small fraction 
of metals in the clouds driven by the galactic fountains is locked 
into dust grains. Can this amount significantly increase by dust 
accretion during the journey of the cloud? To answer this ques- 
tion we have to compare the dynamical timescale of our clouds 
with the typical dust accretion timescale. 

Following the work of Calura et al. (2008) let Xd be the abun- 
dance by mass of the dust and cr the ISM fraction at the time t, 
the quantity a dust — Xd ■ cr represents the normalized mass den- 
sity of the dust at the time t. 

The accretion timescale T accr is given by: 

T accr = T /(l-X d ), (54) 

for to typical values range between 50 Myr and 200 Myr (see 
Dwek 1998). The integration of the equation for the temporal 
evolution of Xd with the initial condition Xd,o — and « 1 
gives: 

t^T 0J [\n(Xd/Xd fi )-Ml)] (55) 

We want to know how much time does it take to increase the 
initially small amount of dust up to a non-negligible value. If we 
calculate for instance the time necessary for accreting 10 times 
the initial dust fraction (e.g. Xd/Xdfl = 10), we get: 

t - 2.3-ro,;. (56) 

This time is greater than the average time of the cloud orbits, 
therefore we can conclude that the depletion of metals into dust 
does not have a important role in the supershell evolution. 

5. Conclusions 

In this paper we have studied the evolution of a supershell pow- 
ered by the energetic feedback of a typical Galactic OB associ- 
ation at various Galactocentric radii. Based on the Kompaneets 
(1960) approximation, we have found analytical solutions of the 
temporal evolution of the supershell and we have established 
criteria for its fragmentation, which can create clouds that are 
thrown out of the disk. Given the self-similar behavior of the 
Kompaneets solutions, the clouds are formed at the same scale 
height (~ 450 pc), irrespective of the number of SNe in the OB 
association or of the Galactocentric radius at which the OB as- 
sociation are. Assuming that the ejecta of the dying stars of the 
OB association instantaneously mix with the supershell, we are 
able to calculate the chemical composition of the clouds and in 
particular their [O/Fe]. 

We have considered four different OB associations (contain- 
ing 10, 50, 100 and 500 SNe, respectively) and three differ- 
ent initial throwing coordinates (4, 8 and 12 Kpc, respectively). 
Once the clouds are formed and can leave the disk, we follow 
their orbits either assuming a purely ballistic model, or introduc- 
ing a viscous force acting between the cloud and the surrounding 
hot halo gas. Our main conclusions can be summarized as fol- 
lows: 



Spitoni et al.: Galactic fountains 



11 



- If the initial metallicity of the OB association is solar, the 
pollution from the dying stars has a negligible effect on the 
chemical composition of the clouds. In particular, the [O/Fe] 
abundance ratio reaches at most ~ 0.025 in a model in which 
the throwing coordinate is Ro = 12 kpc. Only starting from 
very low metallicities (less than 1/100 Z ) it is possible to 
produce a significant enrichment of a-elements. 

- Both in the ballistic and in the viscous interaction models 
the maximum height reached by the clouds is not very large. 
Only for OB associations composed of 500 SNe it is possible 
to throw clouds up to heights larger than 2 kpc above the 
plane of the Galaxy. 

- The range of the cloud orbits is also quite small. The clouds 
are generally directed outwards but the average landing co- 
ordinates differ from the throwing coordinates by ~ 1 kpc at 
most. Only for a throwing coordinate of 12 kpc and an OB 
association made of 500 SNe the ballistic model predicts a 
landing coordinate ~ 2 kpc larger than the throwing one. 

- Models including a viscous interaction between clouds and 
the extra-planar gas predict smaller ranges of the cloud or- 
bits. Indeed, the drag experienced by the cloud brakes it and 
therefore it shortens its journey above the Galactic plane. 

- The HVC Complex C has a mass, velocity and inferred 
height above the plane inconsistent with the results of our 
models. Its oxygen overabundance ([O / 'Fe]=0.12) can be 
reproduced only if we assume a large OB association with 
metallicity 0.1 Z or smaller, therefore its Galactic origin 
cannot be completely ruled out on the basis of its chemical 
composition alone, but the kinematical data suggest a dif- 
ferent formation mechanism. The IVC Arch IV has instead 
velocities and heights above the plane easily reproduced by 
our models, but its [O/Fe] =0.25 is much larger than the one 
of Complex C and it can be explained only by assuming ini- 
tial metallicities smaller than 1/100 Z which are unlikely 
at the present time for the Galactic disk. Therefore, it is un- 
likely that the two studied clouds are originated in a Galactic 
fountain motion. 

Acknowledgements. We thank the referee for the enlightening suggestions. We 
are grateful to A. D'Ercole for having kindly provided a ballistic galactic foun- 
tains code. We also thank F. Calura, G. Cescutti for many useful discussions. 

References 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta 53, 197 
Andrievsky S.M., Bersier D., Kovtyukh V. V. et al 2002a, A&A 384, 140 
Andrievsky S.M., Kovtyukh V. V., Luck R E. et al 2002b, A&A 381, 32 
Andrievsky S. M., Kovtyukh V. V., Luck R.E. et al 2002c, A&A 392, 491 
Andrievsky S. M., Luck R. E., Martin P. et al 2004, A&A 413, 159 
Asplund, M., Grevesse, N., Sauval, A. J. 2005, in "Cosmic Abundances as 

Records of Stellar Evolution and Nucleosynthesis", eds T. G. Barnes III, F. 

N., Bash, ASP Conf. Ser. 336, 25 
Booth, C. M., Theuns, T. 2007, MNRAS, 381, 89 
Barnabe, M., Ciotti, L., Fraternali, F, Sancisi, R., 2006, A&A, 446, 61 
Basu, S., Johnstone, D., Martin, P. G., 1999, ApJ, 516, 843 
Blitz, L., Spergel, D. N„ Teuben, P. J., Hartmann, D., Burton, W. B. 1999, ApJ 

514, 818 

Bluhm, H., de Boer, K, Marggraf, O., Richter, P. 2001, A&A, 367, 299 
Boomsma, R., Oosterloo, T., Fraternali, F, van der Hulst, T., Sancisi, R. 2005, 

in "Extra-planar Gas"Conference, ASP Conf. Series, ed. R. Braun, vol. 331, 

p. 247 

Bregman, J. N. 1980, ApJ, 365, 544 

Calura, F; Pipino, A.; Matteucci, F, 2008, A&A, 479, 669 

Cescutti, G., Matteucci, F, Francois, P., Chiappini, C. 2006, A&A, 462, 943 

Charlton, J. C. & Salpeter, E.E. 1989, ApJ, 346,101 

Chiappini, C, Matteucci, F, Romano, D., 2001, ApJ, 554, 1044 

Collins, J., A., Benjamin, R. A., Rand, R. J., 2002, ApJ, 578, 98 

Cox, D. P., 2005, Annu. Rev. Astro. Astrophys., 43, 337 

de Zeeuw, P. T., Hoogerwerf, R. de Bruijne J.H.J., 1999, ApJ, 1 17, 354 



Draine, B.T., & Salpeter, E.E. 1979, ApJ, 231, 438 
Dwek, E., 1998, ApJ, 501, 643 

Fraternali, F, Oosterloo, T. Sancisi, R. 2004, A&A, 424, 485 
Fraternali, F, Binney, J., 2008, arXiv:0802.0496 accepted for publication in 
MNRAS 

Grevesse, N., & Noels, A. 1993, in Orgin of the Elements, ed. N. Prantzos, E. 

Vangioni-Flam, & M. Casse, (Cambridge: Univ. Press), 15 
Hernquist L.,1990, ApJ, 356, 359 
Houck, J. C, & Bregman, J. N. 1990, ApJ, 352, 506 
Ibata, R. A., Gilmore, G., Irwin, M. J. 1994, Nature 370, 194 
Kompaneets, A. S., 1960, Soviet Phys. Dokl., 5, 46 
Jones, A.P, Tielens, G.G.M., & Hollenbach, D.J. 1996, ApJ, 469, 740 
Lu, L., Sargent, W. L. W., Savage, B. D., Wakker, B. P., Sembach, K. R., 

Oosterloo, T. A. 1998, AJ, 115, 162 
Luck R. E., Gieren W. P., Andrievsky S. M. et al. 2003, A&A, 401, 939 
Mac Low, M. M., McCray, R., 1988. ApJ, 324, 776 
Matteucci, F, & Greggio, L. 1986, A&A, 154, 279 

Matteucci, F, 2001, The Chemical Evolution Of The Galaxy, Kluwer Academic 

Publishers. 
Melioli, 2008, A&a Submitted . 

Miyamoto, M., & Nagai, R., 1975, Publ. Astron. Soc. Japan, 27, 533 

Murphy, E. M., Sembach, K. R., Gibson, B. K, Shull, J. M., Savage, B. D., Roth, 

K. C, Moos, H. W., Green, J. C, York, D. G., Wakker, B. P., 2000, ApJ, 538, 

L35 

Navarro, J. D., Frenk, C. S., White S. D. M., 1996, ApJ, 462, 563 
Oort, J. H. 1970, A&A 7, 381 

Richter, P., Sembach, K. R., Wakker, B. P., Savage, B. D., Tripp, T. M., Murphy, 

E. M., Kalberla, P. M. W., Jenkins, E. B., 2001, ApJ, 559, 318 
Salpeter, E. E., 1955, ApJ, 121, 161 

Sembach, K. R., & Oosterloo T. A. 1998, AJ, 1 15, 162 

Sembach, K. R., Howk, J. C, Savage, B. D., Shull, J. M. 2001, AJ, 121, 992 

Shapiro, P. R., & Field, G. B. 1976, ApJ 205, 762 

Tenorio-Tagle, G, Bodenheimer, P., Rozyczka, M., 1987, A&A, 182, 120 
Tenorio-Tagle, G, Silich, S. A., Kunth, D., Terlevich, E., Terlevich, R., 1999, 

Mon. Not. R. Astron. Soc, 309, 332 
Wakker, B. P., York, D. G., Howk, C, Barentine, J.C., Wilheim, R., Peletier, R. 

F, van Woerden, H., Beers, T. C, Ivezic, Z., Richter, P., Schwarz, U. J., 2007, 
ApJ, submitted 

Wakker, B. P., et al. 1999, Nature, 402, 388 
Woosley, S. E., Weaver, T. A., 1995, ApJ, 101, 181 



6. Online Material 

In the following tables we summarize all our results concerning 
the abundances of oxygen and iron calculated with our model 
for the initial throwing radial coordinates Rq = 4, 8, 12 kpc and 
for metallicities: Z , 0.1 x Z , 0.01 x Z , 10~ 4 x Z . 
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Table 7. Z = Z 8 kpc 



Spitoni et al.: Galactic fountains 



Number 
of SNe 


M,Fe 56 


M, 0ie 


M, ej 


MshellFe 56 


M S hellO lf> 


M,„, 


X*Fe 56 


Xo 16 


[O/Fe] 


10 


0.95 


10.28 


1.16 x 10 2 


678.40 


5.56 x 10 3 


57.91 x 10 4 


1.17 x 10- 3 


9.60 x 10- 3 


1.97 x 10- 4 


50 


3.95 


49.89 


4.77 x 10 2 


681.40 


5.60 x 10 3 


57.94 x 10 4 


1.17 x 10- 3 


9.67 x 10- 3 


1.36 x 10~ 3 


100 


6.99 


97.55 


8.52 x 10 2 


684.44 


5.69 x 10 3 


57.98 x 10 4 


1.18 x 10~ 3 


9.75 x 10- 3 


3.11 x 10~ 3 


500 


24.69 


437.46 


3.18 x 10 3 


702.14 


5.99 x 10 3 


58.22 x 10 4 


1.21 x 10- 3 


1.03 x 10- 2 


1.74 x 10~ 2 



Table 8.Z = 0.1 xZ Q 8kpc 



Number 
of SNe 




M, 0l6 


M, ej 


MshellFe 56 




M, ot 


X*Fe 5e 


X*o 16 




[01 Fe] 


10 


0.95 


9.33 


1.14 x 10 2 


68.70 


5.64 x 10 2 


57.91 x 10 4 


1.19 X 10~ 4 


9.75 x 10" 


-4 


1.17 x 10~ 3 


50 


3.95 


45.48 


4.81 x 10 2 


71.70 


6.01 x 10 2 


57.95 x 10 4 


1.24 x 10~ 4 


1.04 x 10" 


-3 


9.55 x 10- 3 


100 


6.99 


88.97 


8.59 x 10 2 


75.57 


6.44 x 10 2 


57.99 x 10 4 


1.30 x 10- 4 


1.11 x 10" 


-3 


1.71 x 10- 2 


500 


24.69 


401.60 


3.21 x 10 3 


92.43 


9.57 x 10 2 


58.22 x 10 4 


1.59 x 10- 4 


1.64 x 10" 


-3 


1.01 x 10-' 



Table 9. Z = 0.01 x Z 8 kpc 



Number M, Fe56 M, 0u M, ej M shettFess M sheUOl6 M to , X, Fe56 X, 0l6 [O/Fe] 

of SNe 



10 


0.95 


8.96 


1.14 x 10 2 


50 


3.95 


43.67 


4.82 x 10 2 


100 


6.99 


85.72 


8.61 x 10 2 


500 


24.69 


394.21 


3.23 x 10 3 



7.72 64.50 57.91 x 10 4 

10.73 99.20 57.95 x 10 4 

13.77 1.41 x 10 2 57.99 x 10 4 

31.46 4.50 x 10 2 58.22 x 10 4 



1.33 xlO- 5 1.11 x 10 4 8.29 xlO- 3 

1.85 xlO- 5 1.71 xlO" 4 5.23 xlO" 2 

2.37 x 10- 5 2.43 x 10~ 4 9.74 x 10~ 2 

5.40 x 10~ 5 7.72 x 10^ 4 2.41 x 10"' 



Table 10. Z= l(T 4 xZ 8kpc 



Spitoni et al.: Galactic fountains 
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Number M, Fese M» 0l6 M, ej M S h e iiFe 56 M she uo ie M tol X, Fe56 X, 0l6 10/ Fe] 
of SNe 

10 0.95 8.55 1.14xl0 2 1.02 9.10 57.91 x 10 4 2.09 xlO" 6 1.57 xlO" 5 3.68 x 10~ 2 

50 3.95 41.54 4.81 x 10 2 4.02 42.10 57.95 x 10 4 7.26 x lO" 6 7.26 x 10~ 5 1.06x10-' 

100 6.99 81.58 8.60 x 10 2 7.06 82.14 57.99 x 10 4 1.22 xlO" 5 1.42 xlO" 4 1.52x10-' 

500 24.69 377.66 3.23 x 10 3 24.75 3.78 x 10 2 5 8.22 x 10 4 4.25 x 10~ 5 6.52 x 10~ 4 2.70 x 10~' 



Table U.Z = Z 12kpc 

Number M, Fe}f> M, 0l6 M, ej M she , lFe56 M shd ,o {6 M m X, Fe56 X, 0i6 [O/Fe] 

of SNe 

10 0.92 1 0.23 1.13 x 10 2 456.70 3.75 x 10 3 3 8.97 x 10 4 1.17 x 10~ 3 9.61 x 10~ 3 3.13xl0~ 4 

50 3.74 49.28 4.47 x 10 2 459.52 3.78 x 10 3 39.00 x 10 4 1.18 x 10" 3 9.70 x 10" 3 2.14 xlO" 3 

100 6.35 95.73 7.95 x 10 2 462.13 3.83 x 10 3 39.03 x 10 4 1.18 x 10~ 3 9.81 x 10~ 3 4.98 x 10~ 3 

500 23.27 421.95 2.95 x 10 3 479.05 4.16 x 10 3 39.25 x 10 4 1.22 xlO" 3 1.06 xlO~ 2 2.48 x 10" 2 



Table 12. Z = 0.1 xZ 12kpc 

Number M, Fe}f> M, 0l6 M, ej M shellFe56 M she iio l6 M m X, Fe56 X, 0i6 [O/Fe] 
of SNe 

10 0.92 9.29 l.llxlO 2 46.50 3.83 x 10 2 38.97 x 10 4 1.19xl0- 4 9.82xl0~ 4 1.98xlO~ 3 

50 3.74 45.00 4.51 x 10 2 49.32 4.18 x 10 2 39.00 x 10 4 1.26 xl0~ 4 1.07 xlO~ 3 1.51 x 10~ 2 

100 6.35 87.27 8.02 x 10 2 51.92 4.61 x 10 2 39.03 x 10 4 1.33 x 10~ 4 1.18 x 10~ 3 3.46 xlO~ 2 

500 23.27 3 88.11 2.97 x 10 3 68.85 7.62xl0 2 39.25 x 10 4 1.75 x 10 -4 1.94xlO~ 3 1.30X10" 1 
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Table 13. Z = 0.01 x Z 12 kpc 



Spitoni et al.: Galactic fountains 



Number M, Fe5( 
ofSNe 



M, 0u 



M„ 



MshellFeK MshellO,. 



X*Fe 56 



[O/Fe] 



10 
50 
100 
500 



0.92 8.93 1.11 x 10 2 

3.74 43.24 4.52 x 10 2 

6.35 84.64 8.06 x 10 2 

23.27 381.99 3.00 x 10 3 



5.48 46.29 38.97 x 10 4 1.41 x 10" 5 1.19 xlO" 4 1.31 x 10" 2 

8.30 80.60 39.00 x 10 4 2.13 x 10~ 5 2.07 x 10~ 4 7.36 x 10~ 2 

10.90 121.82 39.04 x 10 4 2.79 x lO" 5 3.12 xlO" 4 1.35x10-' 

27.83 419.35 39.25 xl0 4 7.09 x 10~ 5 1.07 xlO" 3 2.64 x 10"' 



Table 14. Z = 10~ 4 x Z 12 kpc 



Number 
of SNe 


M,Fe 5t , 


M t0l6 


M, ej 


MshellFe 5f , 


MshellO l6 


M lo , 


X,Fe 56 


X,o i6 


[O/Fe] 


10 


0.92 


8.51 


1.11 x 10 2 


0.97 


8.88 


38.97 x 10 4 


2.49 x 10- 6 


2.28 x 10- 5 


4.80 x 10" 2 


50 


3.74 


41.13 


4.52 x 10 2 


3.79 


41.5 


39.00 x 10 4 


9.72 x 10~ 6 


1.06 x 10^ 4 


1.26 x 10-' 


100 


6.35 


80.47 


8.05 x 10 2 


6.39 


80.84 


39.04 x 10 4 


1.64 x 10- 5 


2.07 x 10- 4 


1.88 x 10-' 


500 


23.27 


367.43 


2.99 x 10 3 


23.31 


367.81 


39.25 x 10 4 


5.94 x 10- 5 


9.37 x 10- 4 


2.84 x 10-' 



Table 15. Z = Z 4 kpc 



Number 
of SNe 


M, Fe56 


M, 0l6 


M vj 


MshellFe 56 


MshellO i6 


M to , 


X*Fe 56 




X,o 16 


[O/Fe] 




10 


0.97 


10.32 


1.19 x 10 2 


1179.14 


9.67 x 10 3 


100.71 x 10 4 


1.17 X 10" 


-3 


9.60 x 10- 3 


1.04 x 10" 


-4 


50 


4.26 


50.57 


5.20 x 10 2 


1182.42 


9.71 x 10 3 


100.75 x 10 4 


1.17 x 10- 


■3 


9.63 x 10- 3 


7.02 x 10" 


-4 


100 


7.79 


99.37 


9.33 x 10 2 


1185.95 


9.76 x 10 3 


100.79 x 10 4 


1.18 x 10" 


-3 


9.69 x 10~ 3 


1.58 x 10" 


-3 


500 


27.36 


457.84 


3.53 x 10 3 


1205.52 


10.11 x 10 3 


101.05 x 10 4 


1.19 x 10" 


-3 


1.00 x 10- 2 


1.01 x 10" 


-2 



Table 16. Z = 0.1 x Z Q 4 kpc 



Spitoni et al.: Galactic fountains 
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Number 
of SNe 


M tFe56 


M, 0i6 




MshellFe 56 


MshellO if> 


M to , 


X*Fe 56 


X,o l6 


10/Fe] 


10 


0.97 


9.36 


1.16 x 10 2 


118.79 


9.75 x 10 2 


100.71 x 10 4 


1.18 x 10- 4 


9.68 x 10- 4 


6.17 x 10- 4 


50 


4.26 


45.99 


5.19 x 10 2 


122.08 


10.11 x 10 2 


100.75 x 10 4 


1.21 x 10- 4 


1.00 x 10-' 


4.77 x 10~ 3 


100 


7.79 


90.67 


9.41 x 10 2 


125.61 


10.56 x 10 2 


100.79 x 10 4 


1.25 x 10~ 4 


1.05 x 10~ 3 


1.11 x 10~ 2 


500 


27.36 


418.89 


3.56 x 10 3 


145.17 


13.84 x 10 2 


101.05 x 10 4 


1.44 x 10- 4 


1.37 x 10- 3 


6.58 x 10- 2 



Table 17. Z = 0.01 x Z 4 kpc 



Number 
of SNe 


M, Fe56 


M, 0l6 


M, ej 


MshellFe 56 


Mshell0 16 


Mtot 


X*Fe 5t> 


X,o if> 


lO/Fe] 


10 


0.97 


9.00 


1.16 x 10 2 


12.76 


1.05 x 10 2 


100.71 x 10 4 


1.27 x 10~ 5 


1.05 x 10~ 4 


4.06 x 10- 3 


50 


4.26 


44.18 


5.20 x 10 2 


16.04 


1.41 x 10 2 


100.75 x 10 4 


1.59 x 10~ 5 


1.40 x 10^ 4 


2.96 x 10~ 2 


100 


7.79 


87.06 


9.43 x 10 2 


19.57 


1.84 x 10 2 


100.79 x 10 4 


1.94 x 10~ 5 


1.82 x 10~ 4 


5.87 x 10- 2 


500 


27.36 


408.93 


3.58 x 10 3 


39.14 


5.05 x 10 2 


101.06 x 10 4 


3.87 x 10- 5 


5.00 x 10- 4 


1.97 x 10-' 



Table 18. Z = 10~ 4 x Z 4 kpc 



Number 
of SNe 


M, Fes6 


M, 0l6 


M, ej 


M S hellFe 56 


MshellO l6 


M to , 


XtFest, 


X*o l6 


lO/Fe] 


10 


0.97 


8.58 


1.16 x 10 2 


1.09 


9.55 


100.71 x 10 4 


1.08 x 10- 6 


9.48 x 10-" 


2.89 x 10- 2 


50 


4.26 


42.06 


5.20 x 10 2 


4.38 


43.03 


100.75 x 10 4 


4.35 x 10-" 


4.27 x 10- 5 


7.86 x 10- 2 


100 


7.79 


82.80 


9.42 x 10 2 


7.91 


83.77 


100.79 x 10 4 


7.85 x 10-" 


8.31 x 10~ 5 


1.11 x 10-' 


500 


27.36 


390.33 


3.58 x 10 3 


27.47 


391.30 


101.05 x 10 4 


2.71 x 10~ 5 


3.87 x 10~ 4 


2.40 x 10-' 



